Effects of spin density wave quantization on the electrical transport in epitaxial Cr 

thin films 
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We present measurements of the electrical resistivity, p, in epitaxial Cr films of different thicknesses 
grown on MgO (100) substrates, as a function of temperature, T. The p{T) curves display hysteretic 
behavior in certain temperature range, which is film thickness dependent. The hysteresis are related 
to the confinement of quantized incommensurate spin density waves (ISDW) in the film thickness. 
Our important finding is to experimentally show that the temperature Tmid where the ISDW changes 
from A'^ to A'' + 1 nodes decreases as the film thickness increases. Identifying Tmid with a first order 
transition between ISDW states with TV and A'^ + 1 nodes, and using a Landau approach to the 
free energy of the ISDW together with Monte Carlo simulations, we show that the system at high 
temperatures explores all available modes for the ISDW, freezing out in one particular mode at a 
transition temperature that indeed decreases with film thickness, L. The detailed dependence of 
Tmid{L) seems to depend rather strongly on the boundary conditions at the Cr film interfaces. 

PACS numbers: 75.30.Fv, 73.50.-h, 75.70.Ak 



I. INTRODUCTION 

Besides the superconducting state, itinerant antifer- 
romagnetism is another macroscopic manifestation of a 
collective quantum state of the charged electrons in a 
metal-. Bulk Cr is the paradigmatic and unique simple 
metal showing this electronic ground state. In single crys- 
talline samples, an incommensurate spin density wave 
(ISDW) forms below the Neel temperature Tat = 311 K, 
characterized by wave vectors Q± = (27r/acr)(l±'^) along 
the (100) directions^. Here acr- = 2.88 A is the lattice pa- 
rameter of the Cr bcc unit cell and 6 = acr/-^iSDW is 
the deviation from commensurability, with Ajsdw the 
wavelength of the ISDW. This deviation is temperature 
dependent)^ leading to a variation of Ajsdw from 60 A 
at 10 K to 78 A at 310 K. Above the spin-flip tempera- 
ture Tgf = 123 K, the spin waves are transverse with the 
spin S perpendicular to Q. Contrary to what happens 
in single crystals in which all three crystallographic di- 
rections are equivalent, in thin Cr films^ the wave vector 
Q orients perpendicular to the film surface as its thick- 
ness decreases^, and the ISDW's are transverse for all 
temperatures^ below T/v. This confinement of the ISDW 
leads to the quantization of the wave vector Qat, with 
A'^ an integer, giving rise to very interesting hysteretic 
phenomena. 

FuUerton et alA measured neutron scattering and elec- 
trical transport in Cr/Cr97.5Mn2.5 (001) superlattices, 
where the CrMn layers, being antiferromagnetic with a 
higher T^, exert some stringent boundary conditions on 
the Cr surfaces. They found that temperature cycling 
produced irreversibility in the derivative of the resistiv- 
ity and in the neutron scattering scans. These features 
were attributed to a change in an odd number of nodes in 
the ISDW inside the Cr layers which is compatible with 
the magnetic coupling at the interfaces with the CrMn 
layers. 



Recently, Kummamuru et al.~ reported on transport 
results in epitaxially grown Cr (100) thin films of dif- 
ferent thicknesses. Here the boundary conditions are 
different from those in the superlattices of Fullerton et 
al. However, they found thermal hysteretic behavior of 
the Hall coefficient and resistivity for thicknesses below 
510 A. They interpreted these results in terms of ISDW 
domains with N nodes (at high temperatures) or N+1 
nodes (at low temperatures) in the Cr film. For thicker 
films (3500 A) this mechanism disappears, and only an ir- 
reversibility in the p(T) curves was seen in a wider tem- 
perature range, and explained in terms of domain- wall 
scattering of electrons. 

If hysteretic phenomena are indeed related to a first 
order phase transition taking place between domains of 
ISDW's with iV and A^-1-1 nodes inside the Cr layers, 
it is natural to ask how the transition temperature de- 
pends on the Cr thickness. Since for thicker films the 
number of available modes for the ISDW increases in 
the whole temperature range where AjsDwiT) changes, 
one would naively expect to find the hysteretic behav- 
ior at higher temperatures in the thicker films. In this 
article we address these issues both experimentally and 
theoretically, and demonstrate that exactly the opposite 
behavior is found, i.e. the hysteretic region and the tran- 
sition temperature move towards lower temperatures as 
the film thickness is increased. We therefore believe that 
this work throws new light on the understanding of the 
behavior of ISDWs in confined geometries. 



II. EXPERIMENT 

Cr thin films were grown epitaxially on MgO (100) sub- 
strates using DC Magnetron Sputtering in a similar way 
as that reported by Kummamuru et al.i^ The films, with 
thicknesses between 270A and llOOA, were characterized 
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FIG. 1. (Color online) Resistivity as a function of temperature FIG. 2. (Color online) Derivatives of the p{T) curves shown 
for a cooling- warming cycle for the 270 A film. The inset in Fig. [T] The inset shows the difference of dp/dT between 
shows the difference (pc ~ Pw)/pw, expressed as a percentage. the cooling and warming cycles. 



by x-ray diffraction showing rocking curves that have an 
angular dispersion at FWHM around the [002] peak of 
0.5°, indicating a very good degree of epitaxy. Results 
from AFM scans show a mean surface roughness of about 
20 A, which is the third part of Ajsdw at low tempera- 
tures. All films were patterned using photolithography 
and chemical etching in a convenient configuration to 
perform four terminal measurements of the resistance. 
It was measured as a function of T between 15 and 340K 
in a commercial cryo-cooler. The error in the resistance 
measurement was 1/10000, almost independent of T. 



III. RESULTS AND DISCUSSION 

In Figll]we show the resistivity of the 270 A film as a 
function of temperature in the temperature region where 
the hysteretic behavior appears in p{T) while making a 
cooling-warming cycle. The inset shows the difference 
between PcooUng and Pwarmmg {Pc~Pw) normalized by p^, 
and expressed as a percentage. Clearly, in a rather small 
temperature window, the resistivity while cooling down 
is higher than the values of p{T) when warming up the 
film by an amount that reaches 0.5% at it's maximum. 
This value is of the same order of magnitude than that 
reported by Kummamuru et ali^ 

We observed that the derivative of the resistivity 
curves has an intrinsic noise in the temperature region 
where the hysteresis is seen, which is not related to the 
measurement method. We will address this issue in a 
forthcoming article. For the purpose of the physics we 
want to analyze here, it is enough to work with the 
smoothed curves of dp/dT as a function of T as shown 
in FigUl corresponding to the data displayed in FigH] 
Clearly, the hysteretic behavior seen in p{T) is refiected 
also in the derivatives of both, the cooling and warm- 
ing cycles. The inset shows the difference between them. 
We define the temperature at which the derivatives are 



equal within the hysteretic zone (marked with an arrow) 
as Tmid{= 235 K for 270 A). This coincides with the tem- 
perature at the maximum of {pc — Pw)/ Pw The tem- 
peratures where the difference between the derivatives 
has its maximum and minimum values, below and above 
Tmid: are marked in the insets of Figs IT] and [5] using dot- 
ted vertical lines. We call these temperatures Tl and 
Th (212 K and 254 K respectively, for 270 A) and define 
AT = Th -Tl. As we wih show below, the AT region 
where the hysteresis occurs is film thickness dependent. 

To appreciate the mentioned behavior, we plot in Figl3] 
Tmid as a function of film thickness L. Note that the 
bars are not error bars but represent the irreversibility 
region width, AT, as defined above. It can be observed 
that there exists a systematic shift of the hysteretic re- 
gion towards lower temperatures as the film thickness is 
increased. In fact, T^id goes down and tends to sat- 
urate as the film thickness becomes larger than 750 A. 
For comparison, we include the data from Kummamuru 
et alj2,(squares). As can be seen, their data are in very 
good agreement with our results. This would be expected 
because both set of films have been grown following the 
same procedure and over the same type of substrate. In 
the data of FuUerton et al.— a similar decreasing trend 
for Tmid may be observed, but a much steeper decrease 
of Tmid with thickness seems to be present, at least for the 
samples reported (63 A and 200A) . Since it is well known 
that the magnetism at Cr surfaces strongly depends on 
conditions at its surface or interface^, this different be- 
havior might be related to the role played by the CrMn 
layers in their Cr/CrMn superlatticesi^i 

Other important features are worth to remark. In the 
inset of Fig. [3]wc plot the difference of the derivatives as 
a function of T for films of increasing thickness. For the 
sake of clarity, the vertical zero of each curve has been 
shifted by an arbitrary offset. The important point here 
is not the absolute value of dp/dT, but the temperature 
Tmid at which dpc/dT (cooling) equals dp^/dT (warm- 
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FIG. 3. (Color online) Transition temperature as a function 
of film thickness. Red squares correspond to data in ref..^ 
The vertical bars are not error bars but the width of the hys- 
teretic region (see text). The dashed lines are guides to the 
eye. Inset: Difference between the derivatives of p{T) curves 
for cooling and warming cycles of films with increasing thick- 
nesses. For clarity, the vertical zero of each curve (horizontal 
lines) has been arbitrarily shifted. The arrows indicate the 
temperature Tmid (see text). 



ing). Note that AT increases, while the amplitude of the 
difference between dpc/dT and dp^jdT, measured at 
and T//, decreases. These two last observations give a 
clear clue on why the irreversible behavior in p(r) is not 
likely to be seen in thicker samples. 

If we identify Tmid as the temperature at which a 
first order transition takes place between a state with 
an ISDW with iV nodes to a new state with iV+I nodes 
when the temperature is reduced, then our data clearly 
show that the transition temperature decreases as the 
film thickness is increased, tending to saturation for films 
thicker than 750 A. In order to explain this behavior of 
Tmid{L), let us look more in detail the quantization of 
the ISDW in thin films, which is important to under- 
stand the hysteresis of the resistivity in cooling-warming 
cycles. 

As a product of different internal interactions, 6 de- 
pends on temperature.'^ It varies quite linearly from, 
(5 = 0.0379 at T = 300K, to 5 = 0.044 at r = 200K. As a 
direct consequence, Kjsdw changes from 76 A to 65.5 A 
in the same temperature range. At much lower temper- 
atures it saturates to a constant value (5c = 0.048 giving 
A/S_DW«60A. 

Since the magnetic instability in the electron system 
depends on the combination cos(Q+ • r) + cos(Q_ • r), 
it is clear that the resulting beating has a fast variation 
modulated by a long wavelength ((5+ — Q-)/2 ~ 
Thus, the ISDW in Cr has a wavelength 
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When the ISDW is confined in Cr films with non- 



FIG. 4. (Color online) Available ISDW modes for different 
film thicknesses in the temperature region of interest. The 
temperature scale is the same for all panels. 



magnetic interfaces, as is the case in our Cr films, we 
expect a configuration with antinodes (or nodes )^ at the 
films' surfaces. For this case, the available modes are sep- 
arated by half a wavelength. Given a film with thickness 
L, the allowed values for Ajv are therefore 2L/N , where 
N is either the number of nodes inside the film (antin- 
odes at the surfaces), or antinodes in the film (nodes at 
the surface). 

In the window of wavelength values found in bulk 
chromium given by 5{T) (60 A to 78 A), there are only a 
few modes available for each film thickness. Their num- 
ber increases as L increases. This is shown in FigUlwhere 
we plot the allowed modes for each film thickness. The 
lines correspond to A(r) as calculated from equation [TJ 
using a linear approximation for 5{T) between 200 K and 
300 K, obtained from the data in Ref.'^. This approxima- 
tion seems reasonable because in the thickness range in 
which our samples lay (270-1100 A) we do not see appre- 
ciable variations of the Neel temperature, indicating that 
there is no strong dependence of 5 with thickness J^ii^ 

As can be seen in FigH in the the film^ of 175 A, the 
transition occurs between N = 4 and N = 5 at the temper- 
ature where the line corresponding to A5 « ^isdw{T)- 
Correspondingly, for the 270 A film the transition be- 
tween N = 7 and N = 8 lies at a lower temperature. From 
this, it seems that the transition temperature would de- 
crease with L. This reasoning, however, cannot be ex- 
tended to the other film thicknesses based on the in- 
formation displayed in Fig|4] alone, because more modes 
than one become available for the ISDW, and there is 
no way to choose a priori at which one the transition 
would occur. In order to understand how the transition 
temperature depends on L, in Section IIVI we propose a 
phenomenological model which takes into account this 
variation of the number of permitted modes when L in- 
creases, and we also incorporate some dynamics. 
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IV. MODEL 

-1 

A. The Free Energy of the ISDW 



When the temperature is decreased below Tj^, the 
wavelength of the ISDW in the Cr film would try to 
follow it's bulk behavior.'^ However, as discussed in the 
previous section: the ISDW is confined in a film of a 
given thickness, the variation in temperature of Kjsdw 
is quantized, and the number of nodes varies discontin- 
uously from TV to ./V + l, in a process that involves a 
first order phase transition. Neutron diffraction^ as well 
as X-ray experiments^ confirm this physical picture in 
Cr/CrMn multilayers and in Cr thin films respectively. 

Let us consider a one-dimensional model of a trans- 
verse ISDW with Q = qz, where z is the direction per- 
pendicular to the film's surface. In this case the order 
parameter takes the form ^ = ^E'oe*'^^. We can then write 
an expression for the free energy (per unit length) de- 
pending on 'So and g, 

+ C\{V -iqo)-^\^ + U SUV face, (2) 

where is the part of the free energy which is inde- 
pendent of the ISDW, T^""= =311K and qo is the wave 
vector that the ISDW would have in the bulk material. 
Following McMillani^ we include the term |(V — i(7o)*P 
which gives a contribution proportional to (g — (70)^^0- 
This takes into account the fact that there exists an en- 
ergy cost in order to deform the ISDW with respect to 
it's bulk form. The key point is that varies in tem- 
perature and these variations should be followed by the 
ISDW in order to keep the energy at a minimum. 

There is, nevertheless, a competition between the en- 
ergy paid to deform the ISDW respect to the bulk, and 
the energy cost for not satisfying the boundary conditions 
imposed by the surface. As already mentioned, the most 
natural way to choose these boundary conditions for Cr 
films that are not in contact with magnetic materials is to 
have nodes or antinodes at the surfaces. For both cases 
An — '^L/N , i.e. q = nN/L — k^. Considering this, 
the surface contribution to the free energy is given by a 
term proportional to {q — fcAr)^^Q. The way to account 
for the different possible available modes {k^} is simply 
by selecting the one that gives the lowest contribution to 
the free energy at a given temperature (see Fig|4]). Then, 
we explicitly have for the free energy, 

F[^o, q]=Fo + a{T - T^""=)4'2 + b-ff, 

+ ciq-qo)^^l + Umm[{q-kN)^^l]. (3) 
{few} 

If the nodes (or antinodes) are pinned at the surface 
(enough large U), then the free energy displays a series 
of minima corresponding to parabolic potentials sitting 
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FIG. 5. Cuts of the two dimensional free energy, F[5'o, q{T)], 
at different temperatures and for a fixed value of 'i'o- a) 
T/Tn = 0.96, b) T/Tn = 0.74, c) T/Tjv = 0.48. Available 
modes correspond to a film of 760 A. 

close to each k^- This is shown in Fig l5] for a fixed am- 
plitude of the order parameter and adequate numeric 
parameters. Note the change in the potential landscape 
with temperature. As T is decreased the barriers be- 
tween modes increase and the free energy changes in such 
a way to give preference to the modes with higher g's, i.e. 
shorter wavelengths and more nodes (antinodes) in the 
ISDW. 

In this simple way we can understand how the ISDW's 
wavelength is jumping from N to N + I nodes and so on 
as the temperature decreases. Given the fact that each 
of these jumps involves a first order transition, hysteresis 
is expected in cooling- warming cycles. 

B. Monte-Carlo Simulation Results 

One of the most important results of this work is the 
description of how the temperature of the maximum hys- 
teresis, Tmid, moves to lower values as the thickness of 
the Cr films is increased. From the scheme shown in 
FiglHit can be inferred that when many ISDW modes are 
available (larger L's), a large barrier should be needed in 
order to freeze the state of the system with the largest 
number of nodes, otherwise the system would wander be- 
tween all available modes. Since the ISDW's amplitude 
^I'o increases with decreasing temperature as well as the 
barriers do, lower temperatures are needed to freeze the 
system in a high fcjv mode. 

To better quantify this, we introduce some dynamics 
by means of a Monte-Carlo simulation of the system. 
Here we have a two dimensional phase space {q, 4*0) and 
an energy given by Eql21 The resulting average value of 
q after several runs is shown in Figl6] for different film 
thickness as the temperature is decreased. At high tem- 
peratures, the system explores all possible fcjv's for each 
L, but as the temperature is lowered, the ISDW chooses 
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FIG. 6. (Color online) Average g as a function of temper- 
ature for different film thicknesses. The fluctuations in q at 
high temperatures freeze to a fixed number of nodes at a tem- 
perature Tmid which depends on the film thickness (arrows). 
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freezing temperature corresponding to each film thick- 
ness, are plotted in Figl?! Clearly, the freezing temper- 
ature follows the same trend as the experimental data 
for T„iid{L), i.e., it goes down as the film thickness is 
increased. The saturation seen in the transition tem- 
perature in our data for Us larger than 750 A naturally 
occurs in our Monte-Carlo simulations if we put, mim- 
icking the experimental results,'^ a constant 5{T) at low 
temperatures. With this simple phenomenological model 
we do not pretend to get a quantitative agreement with 
our experimental data, but it is important to understand 
the physics behind the observed hysteretic phenomena. 
V. CONCLUSIONS 

We have shown that the hysteretic behavior that ap- 
pears in p{T) in epitaxial Cr thin films is a consequence 
of the quantization of a transverse incommensurate spin 
density wave in the direction perpendicular to the film 
surfaces, with a first order phase transition occurring be- 
tween domains with N and -I- 1 nodes (or antinodes) . 
The number TV at which the transition manifests is film 
thickness dependent, being larger for thicker films. An 
analysis of the one dimensional ISDW free energy, esti- 
mated using McMillan's approach for CDW's, together 
with some dynamics introduced via Monte-Carlo simula- 
tions, show that the potential barriers between allowed 
modes increases as the system is cooled down. As a con- 
sequence of this, and for each i, at high temperatures the 
system explores all available modes but remains frozen in 
a high /cat state at low enough temperatures. Therefore, 
this simple phenomenological model explains in a natural 
way the experimental results which show that the tem- 
perature at which the transition between N and iV -I- 1 
modes takes place, decreases as the film thickness is in- 
creased. 



FIG. 7. Freezing temperature from Fig[6]as a function of film 
thickness. The dashed line is a guide to the eye. 



only one fc^v and stays there no matter how low the tem- 
perature goes. This behavior defines a freezing tempera- 
ture that we identify with Tmid- 

The values obtained from the results in Figl6] for the 
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